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ABSTRACT 

Interpretation of the oscillation spectrum of the slowly pulsating B-type star HD21071 
is presented. We show that non-rotating models cannot account for the two highest 
amplitude frequencies and taking into account the effects of rotation is necessary. 
Rotating seismic models are constructed using various chemical compositions, opacity 
data, core overshooting parameters and rotational velocities. There are prospects for 
seismic modelling of SPB stars, even if no asymptotic pattern is observed in their 
oscillation spectra, provided an unambiguous mode identification is doable and the 
effects of rotation are properly included. 
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1 INTRODUCTION 


Main sequence stars of mid to late B spectral types pulsating 
with periods of the order of days have b een dubbe d slowl y 
pulsating B-type stars (SPB stars) by IWaelkeii^ lll99lll . 
Although multiperiodicity and long periods undoubtedly 
pointed to pulsations in high order gravity modes, the exci¬ 
tation mechanism remained unknown at that time. 

Progress has been made possible soon af- 
ter publication of the revised opacity data OPAL 
jlglesias. Rogers fc Wils^ Il992l l with a new maxi¬ 
mum around logT « 5.3 (the so-called metal or Z 
opacity bump) caused by numerous absorption lines 
of iron-group elements. Using the new opac ity data, 
DzienAowgkj ^^foskaliyfc PamvatnvkhI (Il993l l and 


Gautschv fc Sai^ i 1993ll showed that high radial order 
gravity modes in SPB stars are driven by the classical 
K-mechanism operating in the metal opacity bump. 

However, studying SPB stars is not an easy task from 
both observational and theoretical points of view. Obser¬ 
vations of SPB stars are challenging because their low fre¬ 
quencies demand long runs of observations. Furthermore, 
extracting periods of the order of a day is often compli¬ 
cated because of an aliasing effect, especially from ground- 
based observations. Therefore until the Hipparcos mission 
only a few S PB stars were known. Thanks to the Hipparcos 
photometry. I Waelkens et al.l (Il998ll increased the number of 
SPB stars from 11 to about 100. Since the pulsation periods 
are often of the same order as the rotation periods, the next 
obstacle is to distinguish them from each other. In many 
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cases this task is impossible without detailed spectroscopic 
analysis. 

Theoretical interpretation of the SPB pulsations is 
complicated because of three reasons: mode identification, 
very dense theoretical oscillation spectra and rotation. 
Firstly, because of the lack of clear structures in the 
observed oscillation spectra, the only way to determine 
the spherical harmonic degree, i, and the azimuthal or¬ 
der, m, is the use of the information contained in the 


light 

and line profile variations (e. e. Balona & Stobid 

197 c 

: IWatsonI Il988l: ICugier. Dziembowski & Pamvatnvkh 

1994 

: Campos & Smith! 1980a b: Balona! 1986a b: 

Cugier & Daszviiskal I 2 OOII: ICies & Kullavaniiaval Il988: 

Kennellv & Walked Il996l: iTelting & SchriiversI 1199711. 


An exception are the most recent observational re- 
sults from the CoRo T and Kepler space missio n s (e.g ., 
iDegroote et al.l l2010l : iPapics et al.l |2012| . I20l4 1201^ 1 . 
However, we need to emphasize that in very dense 
oscillation spectra such structures can be accidental. 
An example is the SPB star HD 50230. IDegroote et al.l 
(I 2 OI 0 II claim that in the oscillation spectrum of the star 
there is present a sequence equally-spaced in period 
associated with asymptotic properties while our s tudies 
dSzewczuk. Daszvriska-Daszkiewicz fc Dziembowskil l2014h 
contradict this. Secondly, because of very dense theoretical 
oscillation spectra of SPB stellar models, assignment of 
the radial order, n, to the individual observed peaks is 
usually very ambiguous. The last problem concerns the 
effects of rotation, which influence, both, the equilibrium 
model and pulsational properties. In the case of SPB stars, 
even at slow rotation, pulsational frequencies can be of 
the order of the rotational frequency and the effects of the 
Coriolis force cannot be neglected. On the other hand, the 
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advantage is that the influence of centrifugal deformation 
on l ow frequency gravity modes is small ilBallot et alJ 
l 2012 f) . Effects of rotation on low frequency g modes was 
studied in the fra mework of the so-called tradition al 
approximation (e. g., Lee_^_S^^ Townsendl 2003a R 


2005: lDaszvhska-Daszkiewicz. Dziembowski fc Pamvatnyklrl 


20070 or u sing the truncated expansio n for the eigenfunc¬ 
tions (e. g., iLee fc S^ll989l : lLeeH200l1 ~) . 

Attempts to match the observed frequencies 
to the theoretical ones have been un dert ake n by 
Wajczak. Szewczuk fc Daszvhska-Daszkiewicj ll2012l . 

2013h . who found seismic models which fit two observed 


frequencies with well identified degrees, I, in the two SPB 
stars HD 74560 and HD 182255. However, they assumed 
that the observed frequencies are the axisymmetric modes, 
i. e., m = 0 , and applied the zero-rotation approximation in 
pulsational calculations. 

Without a doubt, the best example of seismic modelling 
of the SPB stars so far is that of KIC 10526294. Using the 
Kepler time-series photometry, IPapics et al.l (l2014l l found 
19 frequency peaks quasi-equally spaced in period, most 
of which were split into triplets. These authors interpreted 
the 19 central peaks according to the asymptotic theory as 
dipole g modes with consecutive high radial orders, n. How¬ 
ever, KIC 10526294 is unique among the SPB stars because 
of its extremely slow rotation. The rotational period, as de¬ 
duced from the rotationally split modes, is equal to 188 days, 
which justified neglecting the effects of rotation in seismic 
modelling. _ _ 

Furthermore, IPapics et al.l ll2015h found 36 frequency 
peaks quasi-equally spaced in period being probably a man¬ 
ifestation of the asymptotic properties in the SPB star 
KIC 776080. The star seems to be even more attractive than 
KIC 10526294 not only because of the longer equally spaced 
modes series but also because of the higher rotation velocity 
of at least 62 km s“^. Detailed seismic modelling is still to 
be done. 

The goal of this paper is to present results of seis¬ 
mic modelling of the SPB star HD 21071. The star has a 
few frequency peaks determined from ground-based pho- 
tometry with a well identified angular numbers [I, m) 


dSzewczuk fc Daszvhska-Daszkiewicj|2015bt l. The effects of 
rotation are included via the traditional approximation. In 
the next section, we present the star and results of earlier 
studies. In Section 3, the oscillation spectrum of the star 
is interpreted. Section 4 is devoted to seismic modelling. In 
the last section we summarize the results and discuss the 
prospect for future studies. 



Figure 1. The position of HD 21071 (black diamond) in the HR 
diagram with la error box (solid line rectangle) and 3a error box 
(dash-dotted line rectangle). Also shown are evolutionary tracks 
passing through the centre and edges of the error boxes. Evolu¬ 
tionary calculations were performed using the Warsaw-New Jer¬ 
sey evolutionary code, OP opacities, AGSS09 chemical mixture, 
the initial hydrogen abundance Xq = 0.7, metallicity Z = 0.0082 
and without overshooting from the convective core. The rotational 
velocity at the ZAMS was Vrot = 50 km s“^. 


llFreire Ferrero et alll20ld). logT.fr = 4.157 JZorec fc Rover] 
I2OI2I) . logT^ff = 4.130 i Lefever et al.|[ 20 lR) . In this paper 
we ad opted logTefr = 4.164 ± 0.007 obtained bv iNiemczural 
(l2003t) from the lUE ultraviolet spectra. This value is ap¬ 
proximately equal to the mean of the values obtained by 
the above mentioned authors. The metallicity , as deter- 
mined from the lAU spectra, is 0.00821 q'qq 32 I Niemczural 
I2OO3I) . To put the star on the HR diagram presented in 
Fig-B we used logL/Err^ = 2.44 4 ± 0.076 derived by 
ISzewczuk fc Daszvhska-Daszkiewiczl ll2015bl ). 

In Fig.m there are also shown the evolutionary t racks 
calculated with the OP opacity tables llSeatorJ |2005|) and 
the latest heavy element mixture of lAsplund et alT ~( 20091 ) 
(hereafter AGSS09). We assumed the initial hydrogen abun¬ 
dance Xo = 0.7, the metallicity Z — 0.0082 and the 
equatorial rotational velocities on the ZAMS equal to the 
projected value, Vrot sin i = 50 km s“^, obtained by 
lAbt. Levato fc Grossol (l2002l ). Overshooting from convective 
core was not taken into account. The Warsaw-New Jersey 
code (e. g. IPamvatnvkh et ^Il998l ) is used in the evolution¬ 
ary calculations throughout the paper. The estimated evolu¬ 
tionary mass for HD 21071 is M = 3.69 Mq and the ranges 
of masses for the Icr and Sa error box are (3.55, 3.84) and 
(3.28, 4.16), respectively. 


2 THE SPB STAR HD 21071 


HD 21071 (HR 1029, HIP 1598 8, V576Per ) is a star of 
brightness V — 6.1 mag (iRee dl I2005I) and spectral 

type B7V llMorgan. Hiltner fc Garrisonlll97ll : iMoolev et al.l 
I2OI3I ). Slightly different spectral clas sifications can be als o 
found in the literature, e. g., B 6 V (I A bt fc Hunter! 1 19621) . 
B5IV llLeshlll9^ and B4/B3 jFreire Ferrero et all 2012h . 

There are many determinations of the effective tem¬ 
perature from photometric calibrat ions as well as spec - 
tral fitting, e. g. , log Tetr = 4.169 llSaffe fc Levatd l2014h . 


logTeff = 4.149 dSilai fc LandstreCT 12014 ). logTeft = 4.212 


3 OSCILLATION FREQUENCIES OF HD 21071 

HD 21 071 was classified as an SPB star bv lWaelkens et 
(Il998t) who found variability with a period of P = 0.84 d in 
the Hipparcos space photometric dat a. In the seve n band 
of the Geneva photometric system, iDe Cat et al.l (l2007l ) 
detected four frequencies: 1/ 1 = 1.18843d~^ , whic h coin¬ 
cidences with this given by IWaelkens et al.l (Il998l ). U 2 = 
1.14934d“^ (also present in the Hipparcos data), 1/3 = 
1.41968d”’^ and 1/4 = 0.95706d“^. In Table[T] we give the 
values of these frequencies and the corresponding amplitudes 
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Figure 2. Oscillation spectrum of HD 21071. In parentheses are 
given I and m, while ’p’ and ’s’ indicate the frequency observed 
in photometry and spectroscopy, respectively. 


in the UBV Geneva filters (columns from 2 to 5). It should 
be mentioned that due to a strong aliasing there is a risk that 
1/3 and i/4 can be mistaken with their aliases. Two dominant 
freque ncies. ui and U 2 , were also detected bv [Andrews et aP 
ll2003l l in the A PT data. Mo reover ui was detected in spec¬ 
troscopic data llDe Catll20dil . 

Identification of the degre e, i, for the freque ncies 
of HD 21071 was performed by IPe Cat et ah] |2003) and 
the res ults are given in the 7th column o f Table lJ Re¬ 
cently, ISzewczuk fc Daszvhska-Daszkiewicd (|2015al ^ in¬ 
cluded the effects of the Coriolis force in the mode iden¬ 
tification. They determined both angular numbers {i, m) 
as well as constrained the rotational velocity (Vrot £ 
(150, 250) [km s“^]). The values of {£,m) are given in 
the last column of Table[T] For completeness, in the 
pe nultimate column we added identification o f I obtained 
by ISzewczuk fc Daszvhska-Daszkiewicd ll2015al f9i with the 
zero-rotation approximation. 

As one can see from Fig. (2] three frequencies HD 21071, 
i. e. vi, V 2 and 4 / 3 , are dipole axisymmetric modes. The fre¬ 
quency vn, can be a dipole axisymmetric or a retrograde 
mode, or a quadrupole retrograde mode. One can also see 
that 1/4, vi and 123 are equally spaced in frequency. If we are 
not in the asymptotic regime of high order acoustic modes, 
which is obviously true in our case, a rotational origin of 
triplets is a most probable explanation for such type of struc¬ 
tures. If we are dealing with a rotationally split triplet, the 
central peak, ui, should be an axisymmetric mode, the left 
side peak, 124, should be a retrograde mode and the right side 
peak, V 3 , a prograde mode. However, this is in contradiction 
with the identification of 1/3 as an axisymmetric mode. 

Another evidence suggesting the non-rotational origin 
of the triplet is its perfect symmetry. The difference between 
the central and left side peaks, ui — U 4 = 0.23137 ± 0, 00009, 
is equal within the errors to the difference between the right 
side and the central peaks, vs — = 0.23125 ± 0,00007. 

Linear pulsational theory predicts an asymmetry and the 
left side peak should be closer to the central one than the 
right side peak (see also Appendix A for more details, only 
in the electronic edition of the journal). 

Because the study of higher order effects of mode cou- 
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4.176 4.17 4.164 4.158 4.152 4.146 

log Tgff 

Figure 3. Theoretical frequencies of dipole modes as a function 
of the effective temperature. The frequencies were calculated in 
the zero-rotation approximation for the model with the param¬ 
eters M = 3.52Mq, logTeff = 4.1571, logL/LQ = 2.277 and 
Z = 0.010. The smallest separation between consecutive radial or¬ 
ders with frequencies above 1 d“^ is indicated by the two headed 
arrow. The observed frequencies are marked by thick horizontal 
lines. 


pling is beyond the scope of this paper, we assumed that the 
triplet in the oscillation spectrum of HD 21071 is accidental. 


4 SEISMIC MODELLING 
4.1 Zero rotation-approximation 

The two highest amplitude frequencies, ui and zz 2 , are sep¬ 
arated from each other by only 0.03909 d“^. Since from 
the mode identification we know they are dipole axisym¬ 
metric modes, such small separation can suggest that they 
are consecutive radial orders. To check this hypothesis, we 
constructed a preliminary grid of pulsational models in the 
framework of the zero-rotation approximation. We calcu¬ 
lated models lying within the 3 ct error box (see Fig.[TJ for 
three values of metallicities, Z — 0.010, 0.015, 0.020, the 
initial hydrogen abundances, Xq = 0.7, OP opacity tables 
and AGSS09 chemical mixture with a step in mass of 0.005 
Mq. No overshooting from convective core has been taken 
into account. It turned out that a minimal separation be¬ 
tween consecutive radial orders of dipole modes in the fre¬ 
quency range appropriate for ui and U 2 (i.e. above 1 d“^) 
is min {u{gn) — iz(g„+i)} = 0.0618 d~^. This minimum the¬ 
oretical separation is almost two times larger than the ob¬ 
served difference between vi and zz 2 and was achieved be¬ 
tween modes gi 3 and gi 4 in the model with M = 3.52 Mq, 
Z = 0.010, logTeff = 4.1571 and logL/LQ = 2.277. In 
Fig-El there is shown the evolution of the frequencies of 
dipole modes from ZAMS to the rightmost edge of the 
3a error box. As one would expect, min {v{gri) — i'{gn +i)} 
(marked by the two headed arrow) coincides with the 
avoided crossing phenomenon. 

As a consequence, to properly interpret the frequencies 
vi and V 2 , one has to include effects of rotation. In the ro¬ 
tating star the eigenvalue is replaced by A which 

for axisymmetric modes is increasing function of Vrot. In¬ 
creasing A imitates the effect of the higher degree, I, and, as 
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Table 1. Frequencies detected in HD 21071 with the corresponding amplitudes of the light (in the Geneva UBV filters) and radial 
velocity variations. The last three columns contain the results of the mode identifications. 


ID 

(d-i) 

A a 

(mmag) 

A ^ 

(mmag) 

A ^ 

(mmag) 

(km s ^) 

Vrot = 0 (km s“l) 

(e, my 

Vrot G (150, 250) (km s"!) 


1.18843 

±0.00001 

34.3 

±0.8 

21.3 

±0.6 

18.5 

±0.6 

3.27 

±0.85 

1 

1 

(T 0) 

^2 

1.14934 

±0.00003 

15.4 

±0.8 

8.8 

±0.7 

7.7 

±0.6 


1; 2; 4 

1; 2 

(T 0) 

V3 

1.41968 

±0.00007 

6.7 

±0.8 

4.1 

±0.6 

3.8 

±0.6 


1; 2 

1; 2 

(1, 0) 


0.95706 

±0.00009 

5.8 

±0.8 

3.3 

±0.6 

3.0 

±0.6 


1; 2; 4 

1; 2 

(1, 0); (1, -1); 

(2, -1); (2, -2) 


° lDe Cat et al.l ll2007t) : toe Caj ll2002h : ^Szewczuk fc Daszvnska-Daszkiewicj ll2015bl l 


a consequence, the oscillation spectrum is more dense for a 
given frequency range. 

4.2 Including rotation 

In this section, we construct seismic models fitting the 
two {ui, V 2 ) and three {vi, V 2 , 1 ^ 3 ) observed frequencies of 
HD 21071. The effects of rotation were included via the tra¬ 
ditional approximation in which the effects of the Coriolis 
force are taken into account whereas centrifugal distortion 
is neglected. Rigid rotation is assumed. 

The method of mode identi fication developed by 
iDaszvhska-Daszkiewicz et al.l (l2015l ) allows, besides deter- 
mination of i and m, to constrain t he rota tional velocity. 
ISzewczuk fc Daszvnska-Daszkiewicj (l2015bl) showed that 
the rotation rate of HD 21071 is in the range Vrot G 
(150, 250) km s“^. We confine our searching of seismic mod¬ 
els to this range of the rotational velocity. Due to a high 
computational cost, it was done in the two steps. 

Firstly, we constructed a preliminary grid of models 
for Xo ~ 0.7, with the OP opacity data, AGSS09 chemi¬ 
cal mixture and without overshooting from the convective 
core (Oov = 0.0). Models were calculated inside the Sa error 
box (see Fig.[T]) for different masses, metallicities and rota¬ 
tion velocities with a step AM = 0.01 Mq, AZ = 0.00025 
and AVrot = 10 km s“^. Then, to find the approximate pa¬ 
rameters of models which fit the two observed frequencies, 
the theoretical frequencies were interpolated by means of 
multi-dimensional linear interpolation. 

In the second step, around the approximate model pa¬ 
rameters for which theoretical frequencies fit the observed 
ones, we computed a denser grid of models with steps AM = 
0.002 Mq and AZ — 0.00005. A step in effective tempera¬ 
ture in both grids was dynamically determined by the evo¬ 
lutionary code and equal approximately to ATefj ~ 0.0005. 

We will call the grid of models with the input (Aq = 0.7, 
OP, AGSS09, aov = 0.0) as Gl. 

In Fig.[4l we put seismic models fitting vi and 1/2 in 
the HR diagram for the seven values of Vrot- The colours 
are assigned to the values of metallicity, Z. In this figure 
all modes, stable and unstable, were plotted. The use of 



I09 Teff 

0.005 0.01 0.015 0.02 0.025 

z 

Figure 5. Models fitting ui and 1/2 on the diagram log Teff vs. 
M. There are shown models from the Gl grid with Vrot = 150 
km and in which the modes gi4 and gis fit the two observed 
frequencies of HD 21071. 


the instability condition significantly reduces the number 
of seismic models. Our seismic modelling showed that for 
Vrot ^ 220 km s“^ there are very few models with unstable 
modes which reproduce and V 2 . Moreover, these models 
have rather high metallicity. For example, in the case of 
Vrot = 250 km s“^ unstable modes reproducing vi and V 2 
were found in models with Z > 0.02. On the other hand, 
in the case of the lowest considered rotational velocity, ie., 
Vrot = 150 km s“^ unstable modes reproducing vi and V 2 
were found in models with metallicity as low as Z « 0.0054. 

All our seismic models of HD21071 can be divided into 
two families. The first one (hereafter FI), associated with 
the radial orders from n = 13 to n = 19, has on aver¬ 
age higher metallicity and is more luminous. The second 
one (hereafter F2), associated with the radial orders from 
n = 26 to n = 31, has on average lower metallicity and is 
less luminous. In both cases the models which fit zzi and V 2 
become more luminous with increasing rotation velocity. For 
the lowest rates of rotation, i. e., Vrot = 150 and 160 km s“^, 
we found only models with the FI solution within 3cr error 
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4.180 4.170 4.160 4.150 4.180 4.170 4.160 4.150 


log Tgff 



0.005 0.01 0.015 0.02 0.025 

z 

Figure 4. The HR diagrams with the models from the G1 grid which fit vi and 1/2, marked, for seven values of the rotational velocity, 
Vrot. Both, stable and unstable modes were plotted. In each panel, the inner frame corresponds the Icr error box and the entire frame 
to the 3 (t error box. In the bottom right panel there are shown the evolutionary tracks for the three values of metallicity. 
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log Teff 
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Figure 6. The same as in Fig.[4|but models fitting the three frequencies vi, U 2 and 1/3 of HD 21071 are shown. 
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Table 2. The radial orders, n, of dipole modes which fit the three frequencies, ii\, V 2 and 1 / 3 , of HD21071 in the models Gl, G2, G3 
and G4. In square brackets are stable modes. The first column contains the value of the rotational velocity. 



Gl 

G2 

G3 

G4 

Vrot 

OP 

OP 

OP 

OPAL 

(kms“^) 

Xq = 0.70 

Xo = 0.75 

Xo = 0.70 

Xq = 0.70 


ctov = 0.0 

Qov = 0.0 

Oov = 0.2 

ctov = 0.0 

150 

[916: 917, 912] 

[916, 917, 912] 

914, 915, 910 

[916, 917, 912] 


914, 915, 910 



[915, 916, 911] 





914, 915, 910 

160 

913, 914, 99 

913, 914, 99 

914, 915, 910 

- 




913, 914, 99 


170 

914, 915, 910 

914, 915, 910 

[914, 915, 910] 

- 

180 

914, 915, 910 

[915, 916, 911] 

- 

[929, 931, 920] 



914, 915, 910 


914, 915, 910 

190 

- 

- 

- 

[914, 915, 910] 

200 

- 

- 

- 

- 

210 

[929, 931, 920] 

- 

- 

- 

220 

[929, 931, 920] 

[929, 931, 920] 

- 

- 

230 

- 

- 

- 

[929, 931, 920] 

240 

- 

- 

- 

- 

250 

- 

- 

- 

[928, 930, 919] 


box. The seismic models from the F2 solution appear close to 
the lowest edge of the 3cr error box (logL/LQ « 2.22) from 
Vrot = 170 km s~^. With increasing Vrot, we obtained more 
and more models from F2 and less from FI in the space of 
parameters we considered. For the highest rotational veloc¬ 
ity, Vrot = 250 km s“^, we found only seismic models from 
F2. In both cases, there is a clear trend: the lower metallicity 
the higher luminosity. 

The seismic models of HD 21071 lie along sloped lines 
in the HR diagram. For one pair of the radial orders we 
usually have a few nearly parallel sequences of such models 
differing in metallicity. A similar structure occurs in the di¬ 
agram M vs. logTeff. This is because the frequencies of the 
g modes in the range of parameters and frequency we con¬ 
sider, in general, decrease with decreasing effective tempera¬ 
ture and increasing mass. The dependence on Z is varied, for 
the higher radial orders, n > 10 , the pulsational frequency is 
an increasing function of Z, whereas for n < 10 there is an 
oscillatory character of v{Z). In Fig.[Sl for clarity, we show 
the diagram M vs. logTefr for only one pair of the radial 
orders, gi 4 and gis, and one value of the rotational velocity. 

Gaps which sometimes appear along the lines of seismic 
models result from the adopted steps in M, Z, Tefi and 
Vrot, which, although small, can cause the omission of some 
models. 

It Is worth to notice that in the case of FI we have al¬ 
ways consecutive radial orders whereas in the F2 solution 
we have every second radial order. This is because with in¬ 
creasing Vrot the eigenvalue A of the i = 1, m = 0 mode 
increases and as a result of denser oscillation spectrum of 
the high radial order modes is shifted towards higher fre¬ 
quencies. This explains why we have more models from FI 


than F2 for lower rotation rates and vice versa. On the other 
hand coexistence solutions from FI and F2 for fixed rota¬ 
tional velocity (e. g., Vrot = 200 km s“^) is associated with 
considerably different parameters of the models from both 
families of solutions (mainly metallicity and mass) as de¬ 
scribed in previous paragraph. 

An interesting fact is that most unstable modes comes 
from the FI seismic models. The only exception is for 
Vrot = 250 km s“^ where some modes from the F2 solu¬ 
tion are unstable. With increasing Vrot there are less and 
less seismic models with unstable modes fitting the two fre¬ 
quencies ui and U 2 . On the one hand this is due to lower 
metallicity of the F2 models which dominate for the higher 
rotational velocities, on the other, for higher Vrot the insta¬ 
bility domain is shifted towards higher frequencies than the 
observed ones. 

In the next step we selected models which fit the third 
observational frequency, r'a, which also corresponds to a 
dipole axisymmetric mode. Models fitting vi, V 2 and 1/3 are 
shown in Fig.| 6 ] We can see a further reduction in the num¬ 
ber of models. Moreover, we did not find a solution for each 
value of Vrot- The radial orders of modes which fit vi, V 2 
and vz from the Gl grid are given in the second column of 
TableO Unstable solutions exist only for Kot ^ 180 km s“^. 

4.3 Effects of hydrogen abundance, core 
overshooting and opacities 

There are many input parameters, both, from a model and 
microphysics, that may affect pulsational frequencies. Here, 
we examine the influence of the initial hydrogen abundance, 
Ao, overshooting from the convective core (in terms of the 
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parameter Oov) and the opacities. The effect of the core 
overshooting was included ac c ordin g to the formulation of 
iDziembowski fc PamvatnvkhI (l2008ll which takes into ac¬ 
count, both, the distance of the overshooting and partial 
mixing in the overshoot layer. 

To this end we constructed three additional grids of 
models in the same way as explained in Section 4.2. In the 
second grid, G2, we used Xq = 0.75 (comparing to Xq — 0.70 
in Gl), in the third grid, G3, we added the core overshoot- 
ing, Qqv = 0.2, and in the fourth one, G4, we used OPAL 
lllglesias fc Rogerall996l l opacity tables (instead of OP used 
in Gl). 

We found that independently of the adopted grid, all 
trends noted in Gl are also present in G2, G3 and G4. 
Firstly, there are always the two families of solutions. The 
first one, FI, with the radial orders of about n = 15 and the 
second one, F2, with the radial orders of about n — 30. 

With an increased abundance of hydrogen (the G2 
grid), we found less seismic models compared to Gl, in par¬ 
ticular for the F2 solution. 

Seismic models with and without overshooting from the 
convective core (G3 vs. Gl) lie approximately in the same 
places of the HR diagram but the G3 models are confined 
to slightly narrower bands. A similar picture emerges when 
we considered seismic models only with unstable modes. 

Changing the opacity tables (G4 vs. Gl) has very little 
effect on the position of seismic models in the HR diagram. 
The G4 models are only slightly more luminous. The same 
is true if only unstable modes are considered. The exception 
is for seismic models with Vrot = 240 km s“^. In this case, 
unstable modes occurs only in the F2 solution. This is oppo¬ 
site to seismic models calculated with the OP tables. There 
are fewer seismic models with unstable modes obtained with 
the use of the OPAL opacities. This result is not surprising 
as it is well known that for SPB models the computations 
with the OP tables give more unstable mo des than those 
with the OPAL ones fe.g. |Pamvatnvkhlll999l L 

Models which fit the three frequencies vi, 1^2 and vs in 
the grids G2, G3 and G4 are presented in Fig.[3 In Table[^ 
we list the allowed combinations of the radial orders. 

As in the case of Gl, only for the lower rotation velocity, 
Vrot ^ 180 km s“^, we were able to find models which fit 
the three frequencies with the instability condition fulfilled. 
For Vrot 190 km s“^, the number of these seismic models 
is significantly lower and the modes are stable (cf. Table[2]). 

To summarize our results of seismic modelling of the 
star HD 21071, in Table[3] we give the ranges of the stel¬ 
lar parameters of models with unstable modes which fit the 
three frequencies, vi, V 2 and izs, for the four grids. The pa¬ 
rameters are provided for each separate group of models 
with similar properties, shown in Figs. [B] and [7] by ovals. The 
groups are marked as ’grp’ in Tabled In addition, we give 
also the parameters of representative models for each group 
(marked as ’rep’ in Table[31). In the last column we give, for 
the representative seismic models, the goodness of the fit 
defined as 



where v'^ and are observed and theoretical frequency. 


respectively, and an are the observational errors of the fre¬ 
quencies. 

It should be mentioned that masses of seismic mod¬ 
els from Tabled are in the range M G (3.909, 4.661) Mq 
whereas masses of evolutionary models which fall into the 3 ct 
error box are in a slightly wider range M € (3.04, 4.99) Mq. 
The range of masses of evolutionary models given above 
is wider than the one given in Section 2 because here we 
take into account models with all considered metallicities 
{Z € (0.005, 0.025)) and from all grids. 

In Fig.[8l we compare the oscillation spectrum of 
HD 21071 with the theoretical spectra corresponding to the 
representative seismic models from Table[31 In the bottom 
panel of Fig. [HI we show a zoom-in on the frequencies around 
the observed values. As one can see, the model G2g repro¬ 
duces also the frequency within the observational errors 
if it is the mode (. = 1 , m = 0, g22, a possibility consistent 
with our mode identifications. Similarly, the quadrupole ret¬ 
rograde mode £ = 2, m = —2, gas in the model G3i fits the 
frequency V 4 ,. This mode is also allowed by our mode iden¬ 
tifications. 

We can also see that the frequency of the mode I = 
1 , m = 0, g22 is close to V 4 in the models Gld and G4m. 
In a similar distance is also the mode, 1 = 2, m = —2, 
gai, in the model G3k. Given the numerical accuracy of the¬ 
oretical frequencies, which is about Xv = 0.0001 d~^, we 
have to accept also these solutions. The important result is 
that these identifications are valid for the all seismic models 
constructed in this paper. 


5 CONCLUSIONS 

This paper focused on the challenges and prospects for seis¬ 
mic modelling of slowly pulsating B-type stars with the 
effects of rotation taken into account. We used as an ex¬ 
ample the star HD 21071 which pulsates in four frequen¬ 
cies of which the three have a unique mode identification. 
Firstly, we showed that non-rotating models cannot account 
for the two highest amplitude frequencies, vi = 1.18843 and 
U 2 = 1.14934 d“^, and including the effects of rotation on 
high-order g modes is indispensable. Then, having unam¬ 
biguous determination of the two angular numbers {£, m) 
for the three observed frequencies of HD 21071 (zzi, V 2 and 
V 3 = 1.41968 d“^), and constraints on the range of the rota¬ 
tional velocity, Vrot G (150, 250) km s“^, we constructed ro¬ 
tating seismic models which reproduce these three frequen¬ 
cies. We examined the effects of the initial abundance of hy¬ 
drogen, an amount of the core overshooting and the opacity 
data, considering four grids of parameters. Due to the high 
density of the theoretical oscillation spectra, a large number 
of solutions have been obtained. Despite of that, only two 
combinations of the radial orders, n, were allowed: one set 
is around n = 15 and the second one around n = 30. More¬ 
over, the instability condition reduced the number of seismic 
models significantly. Therefore, accurate calculations of the 
opacity data are of the utmost importance. 

Among seismic models fitting the three frequencies, 
we found some that reproduce also the forth frequency, 
1/4 = 0.95706 d“^. In all grids of models only the dipole ax- 
isymmetric modes g22 or quadrupole modes 1 = 2 , m = — 2 , 
g38, g3i have frequencies close to V 4 (given the observa- 
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Figure 7. The same as in Fig. [6] but selected models from the grids G2, G3 and G4, are shown. 
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V4 V2 Vi V3 



Figure 8. The frequencies of dipole axisymmetric modes in the representative models (see Tabled which reproduce three well identified 
frequencies, ui, U 2 and 1^3, of HD 21071. A zoom-in on the frequencies around the observed values is shown in the bottom panel for 
clarity. In this panel the frequency of the quadrupole retrograde mode (£ = 2, m = —2) is added as an inverted triangle. The vertical 
solid lines indicate the observed frequencies and the vertical dashed lines, the observational errors of the frequencies. 
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Table 3. Ranges of the astrophysical parameters for the groups of seismic models (rows indicated by ‘grp’ in the sixth column) and the 
parameters of the selected models for each group (rows indicated by ‘rep’ in the sixth column). All models reproduce the three observed 
frequencies with the unstable dipole axisymmetric modes. Groups of models are defined as models with similar parameters forming 
distinct groups in Figs.|3and[7] In the last column the goodness of fit is given for representative models (cf. Eq.[T}. 


Grid 

opacity 

Ao 


Vrot 

(km s“^) 

ID 

z 

M (Mq) 

logTeff 

log L/Lq 

d 

G1 

OP 

0.70 

0.0 

150 

grp a 

0.0157-0.0160 

4.031 - 4.039 

4.1542 -4.1559 

2.452 - 2.458 

- 

G1 

OP 

0.70 

0.0 

150 

rep a 

0.01575 

4.0357 

4.15528 

2.4554 

0.139 

G1 

OP 

0.70 

0.0 

160 

grp b 

0.0231 - 0.0236 

4.372 - 4.388 

4.1568-4.1597 

2.479 - 2.489 

- 

G1 

OP 

0.70 

0.0 

160 

rep b 

0.02340 

4.3793 

4.15814 

2.4831 

0.331 

G1 

OP 

0.70 

0.0 

170 

grp c 

0.0161 - 0.0164 

4.218 - 4.232 

4.1587-4.1611 

2.541 - 2.549 

- 

G1 

OP 

0.70 

0.0 

170 

rep c 

0.01630 

4.2221 

4.15948 

2.5435 

0.312 

G1 

OP 

0.70 

0.0 

180 

grp d 

0.0216 - 0.0222 

4.351 - 4.364 

4.1547-4.1574 

2.503 - 2.513 

- 

G1 

OP 

0.70 

0.0 

180 

rep d 

0.02210 

4.3521 

4.15492 

2.5042 

0.624 

G2 

OP 

0.75 

0.0 

160 

grp e 

0.0184-0.0188 

4.647-4.661 

4.1686 -4.1712 

2.524 - 2.534 

- 

G2 

OP 

0.75 

0.0 

160 

rep e 

0.01855 

4.6566 

4.17028 

2.5304 

0.230 

G2 

OP 

0.75 

0.0 

170 

grp f 

0.0136 - 0.0139 

4.455 - 4.469 

4.1639 -4.1660 

2.562 - 2.569 

- 

G2 

OP 

0.75 

0.0 

170 

rep f 

0.01375 

4.4643 

4.16493 

2.5657 

0.720 

G2 

OP 

0.75 

0.0 

180 

grp g 

0.0181 - 0.0186 

4.595 - 4.612 

4.1607-4.1636 

2.529 - 2.539 

- 

G2 

OP 

0.75 

0.0 

180 

rep g 

0.01835 

4.6040 

4.16237 

2.5347 

0.003 

G3 

OP 

0.70 

0.2 

150 

grp h 

0.0163 - 0.0165 

3.909 - 3.919 

4.1448 -4.1467 

2.409 - 2.415 

- 

G3 

OP 

0.70 

0.2 

150 

rep h 

0.01640 

3.9139 

4.14568 

2.4114 

0.740 

G3 

OP 

0.70 

0.2 

160 

grp i 

0.0132 - 0.0135 

4.035 - 4.045 

4.1619-4.1639 

2.523 - 2.531 

- 

G3 

OP 

0.70 

0.2 

160 

rep i 

0.01340 

4.0365 

4.16247 

2.5255 

0.497 

G3 

OP 

0.70 

0.2 

160 

grp j 

0.0173 - 0.0176 

4.097-4.106 

4.1571 -4.1587 

2.457 - 2.463 

- 

G3 

OP 

0.70 

0.2 

160 

rep j 

0.01740 

4.0997 

4.15804 

2.4604 

0.214 

G3 

OP 

0.70 

0.2 

160 

grp k 

0.0203 - 0.0213 

4.382 - 4.412 

4.1675 -4.1732 

2.504 - 2.524 

- 

G3 

OP 

0.70 

0.2 

160 

rep k 

0.02115 

4.3860 

4.16833 

2.5072 

0.326 

G4 

OPAL 

0.70 

0.2 

150 

grp 1 

0.0204 - 0.0209 

4.058 - 4.071 

4.1496-4.1523 

2.386 - 2.395 

- 

G4 

OPAL 

0.70 

0.2 

150 

rep 1 

0.02075 

4.0651 

4.15081 

2.3900 

0.413 

G4 

OPAL 

0.70 

0.2 

180 

grp m 

0.0171 - 0.0175 

4.484 - 4.489 

4.1821 -4.1835 

2.609 - 2.614 

- 

G4 

OPAL 

0.70 

0.2 

180 

rep m 

0.01745 

4.4846 

4.18217 

2.6091 

0.433 


tional errors and numerical accuracy). Both modes are in 
agreement with the mode identification. Fitting the forth 
frequency further reduced the number of seismic models. 
However, to derive more stringent constraints on opacities, 
core overshooting etc. more well identified frequencies are 
necessary. Because we tried to constrain seven parameters 
(M, Teff, ^0, Z, K, (Tov, Vrot), we need at least seven well 
identified pulsational frequencies. 

There are many SPB stars with far more rich oscilla¬ 
tion spectra obtained from the space missions, like MOST, 
CoRoT and Kepler, but the problem is their mode identi¬ 
fication. Even if there are some patterns in the oscillation 
spectra, they can be accidental, i.e., composed of modes with 
different pairs {£, m), and cannot be interp reted in term of 
the asymptotic theory dSzewczuk et al.l|2014l ). Therefore, ex¬ 
cept for the excellent example of KIC 10 526294 with 19 con - 
secutive rotationally split dipole modes (IPanics et al.l[20l3) . 
there is no other reliable seismic modelling of an SPB star. 

The aim of this paper was to show that seismic mod¬ 
elling of SPB stars is feasible even if no asymptotic pattern 
is observed in the oscillation spectra and a limited number of 
modes is detected, provided the angular numbers {£, m) are 
determined unambiguously and the effects of rotation are 
properly included. Therefore, a hope is pinned on the two 
colour data from the BRITE mission which, besides more 
frequencies and better accuracy, can be used to identify the 
pulsational modes. 
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APPENDIX A: ASYMMETRY OF THE 
ROTATIONALLY SPLIT MODES IN THE 
FRAMEWORK OF THE TRADITIONAL 
APPROXIMATION 

Linear theory of pulsation predicts asymmetry between ro- 
tationally split components of the given mode, ie., retro¬ 
grade modes should be closer to the axisymmetric ones than 
prograde modes. In Table lXl] we list frequencies of dipole 
modes and the differences between their components for 
the model at the centre of the error box {M = 3.69 Mq, 
logTeff = 4.164, logL/L0 = 2.444, Xo = 0.7, Z = 0.0082, 
R = 2.62 Rq). The rotation splitting was computed with the 
traditional approximation and the rotational velocities were 
chosen to approximately reproduce the value of the observed 
splitting in HD 21071. As one can see, in all cases retrograde 
modes are closer to axisymmetric modes than prograde ones 
and asymmetry is of the order of 0.01 d“^. This property 
is also clearly seen in Fig. lATl where the observed triplet is 
compared with the theoretical spectrum calculated for the 
model at the centre of error box of HD 21071. The observed 
frequencies were shifted in order to align the central peak, 
i^i, with the gi2 (m = 0) mode. 
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Table Al. Theoretical frequencies of dipole modes for several radial orders, n, and the differences between their rotationally split 
components, for the central model of HD 21071 for the three values of the rotational velocity. 


Vrot 

(kms“^) 

n 

u (m = —1) 

(d-i) 

V {m = 0) 
(d-l) 

1/ (m = 1) 

(d-i) 

II 

g 

1 

o' 

II 

u (m = 1) — u (m = 0) 

(d-i) 

60 

13 

0.86855 

1.08271 

1.30635 

0.21416 

0.22364 


12 

0.93075 

1.14612 

1.37041 

0.21537 

0.22429 


11 

1.01204 

1.22871 

1.45373 

0.21667 

0.22502 


10 

1.11148 

1.32946 

1.55497 

0.21798 

0.22551 

62 

13 

0.86667 

1.08742 

1.31828 

0.22075 

0.23086 


12 

0.92855 

1.15066 

1.38221 

0.22211 

0.23155 


11 

1.00942 

1.23300 

1.46530 

0.22358 

0.23230 


10 

1.10854 

1.33333 

1.56619 

0.22479 

0.23286 

65 

13 

0.86408 

1.09459 

1.33623 

0.23051 

0.24164 


12 

0.92547 

1.15755 

1.39986 

0.23208 

0.24231 


11 

1.00584 

1.23941 

1.48263 

0.23357 

0.24322 


10 

1.10431 

1.33940 

1.58344 

0.23509 

0.24404 



Figure Al. The theoretical frequencies (dashed lines) of dipole 
modes calculated with the traditional approximation for the 
model at the centre of the error box of HD 21071 with Vrot = 63 
km s~^. The mode gi 2 is highlighted with a thin solid line. The 
observed triplet (thick solid lines) is shifted by - 0.036 d“^ to 
align its central frequency with the 5-12 {m = 0) mode. 
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